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In this study, the fluctuation-dissipation theory is invoked to shed light on input-output interindus- 
trial relations at a macroscopic level by its application to IIP (indices of industrial production) data 
for Japan. Statistical noise arising from finiteness of the time series data is carefully removed by 
making use of the random matrix theory in an eigenvalue analysis of the correlation matrix; as a re- 
sult, two dominant eigenmodes are detected. Our previous study successfully used these two modes 
to demonstrate the existence of intrinsic business cycles. Here a correlation matrix constructed 
from the two modes describes genuine interindustrial correlations in a statistically meaningful way. 
Further it enables us to quantitatively discuss the relationship between shipments of final demand 
goods and production of intermediate goods in a linear response framework. We also investigate 
distinctive external stimuli for the Japanese economy exerted by the current global economic crisis. 
These stimuli are derived from residuals of moving average fluctuations of the IIP remaining after 
subtracting the long-period components arising from inherent business cycles. The observation re- 
veals that the fluctuation-dissipation theory is applicable to an economic system that is supposed 
to be far from physical equilibrium. 



I. INTRODUCTION 

Both the business cycle and the interindustrial rela- 
tionship are long-standing basic issues in the field of 
macroeconomics, and they have been addressed by a 
number of economists. Recently, we analyzed [1] busi- 
ness cycles in Japan using indices of industrial produc- 
tion (IIP), an economic indicator that measures current 
conditions of production activities throughout the nation 
on a monthly basis. Careful noise elimination enabled 
us to extract business cycles with periods of 40 and 60 
months that were hidden behind complicated stochastic 
behaviors of the indices. 

In this accompanying paper, we focus our attention 
on interindustrial relationship by analyzing IIP data in a 
framework of the linear response theory; the fluctuation- 
dissipation theory plays a vital role in the analysis of IIP 
data. We also discuss the difference between moving av- 
erage fluctuations in the original data and long-period 
components arising from inherent business cycles. The 
residuals may be interpretable as a sign of external stim- 
uli to the economic system. The recent worldwide reces- 
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sion offers us a good opportunity to conduct this study, 
because it delivered an unprecedented shock to the eco- 
nomic system of Japan. 

The interindustrial relations of an economy are con- 
ventionally represented by a matrix in which each col- 
umn lists the monetary value of an industry's inputs and 
each row lists the value of the industry's outputs, in- 
cluding final demand for consumption. Such a matrix, 
called the input-output table, was developed by Leon- 
tief [21 13] . This table thus measures how many goods of 
one industrial sector are used as inputs for production 
of goods by other industrial sectors and also the extent 
to which internal production activities are influenced by 
change in final demand. Leontief 's input-output analysis 
can be regarded as a simplified model of Walras's general 
equilibrium theory [4] to implement real economic data 
for carrying out an empirical analysis of such economic 
interactions. Currently, the basic input-output table is 
constructed every 5 years according to the System of Na- 
tional Accounts (SNA) by the Ministry of Internal Affairs 
and Communications in Japan. 

It should be noted that the input-output table de- 
scribes yearly averaged interindustrial relations. Al- 
though such a poor time resolution of the table may 
be tolerable for budgeting of the government on an an- 
nual basis, various day-to-day issues faced by practition- 
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ers require them to react promptly. We thus need a 
more elaborate methodology that enables investigation of 
the input-output interindustrial relationship with a much 
higher time resolution. 

Econophysics [SHE] is a newly emerging discipline in 
which physical ideas and methodologies are applied for 
understanding a wide variety of complex phenomena in 
economics. We adopt this approach to address the above- 
mentioned issues in macroeconomics. That is, we pay 
maximum attention to real data while drawing any con- 
clusions. Of course it is important to remember that 
real data are possibly contaminated with various kinds of 
noise. The random matrix theory (RMT), combined with 
principal component analysis, has been used successfully 
to extract genuine correlations between different stocks 
hidden behind complicated noisy market behavior [5^16j . 
Recently, dynamical correlations in time series data of 
stock prices have been analyzed by combining Fourier 
analysis with the RMT [T7] . 

In this study, we further develop the noise elimination 
method initiated in previous studies. The null hypothesis 
that has been adopted thus far for extracting true mutual 
correlations corresponds to shuffling time series data in a 
completely random manner. Although we should distin- 
guish between mutual correlations and autocorrelations, 
both these correlations are destroyed at the same time by 
the completely random shuffling. To solve this problem, 
rotational random shuffling of data in the time direc- 
tion is introduced as an alternative null hypothesis; such 
randomization preserves autocorrelations involved in the 
original data. The new null hypothesis thus elucidates 
the concept of noise elimination for mutual correlations. 

Further, we borrow the concept of the fluctuation- 
dissipation theory jl8] from physics to elucidate the in- 
terindustrial relationship and the response of an eco- 
nomic system to external stimuli. The theory establishes 
a direct relationship between the fluctuation properties of 
a system in equilibrium and its linear response properties. 
We assume that the validity of the fluctuation-dissipation 
theory in physical systems is also true for such an exotic 
system as described by the IIP. Very recently, dynam- 
ics of the macroeconomy has been studied in the linear 
response theory by taking an explicit account of hetero- 
geneity of microeconomic agents |19j . 

The present paper is organized as follows. In Sec.|lTj we 
first provide a brief review of the noise elimination from 
the IIP using the RMT. A new null hypothesis based on 
rotational random shuffling is introduced in Sec. |III[ Sec- 
tion |IV| presents construction of a genuine correlation ma- 
trix for the IIP by consideration of only those dominant 
modes that are approved to be statistically meaningful 
by the RMT. We present development of a fluctuation- 
dissipation theory for input-output interindustrial rela- 
tions in Sec. |Vj Then, in Sec. |VI[ we quantitatively 
discuss relationship between shipments of final demand 
goods and production of intermediate goods. In Sec. |VII[ 
we elucidate response of the industrial activities to exter- 
nal stimuli by subtracting long-period components aris- 
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FIG. 1. Input-output relationship in industrial activities of 
economic system as measured by IIP in Japan. The numbers 
in the parentheses denote the classification index g in Tablejl] 



ing from inherent business cycles from moving average 
fluctuations in the original data. Section [VlII| concludes 
this paper. 



II. APPLICATION OF RANDOM MATRIX 
THEORY TO IIP 

In Japan, the IIP are announced monthly by the Min- 
istry of Economy, Trade and Industry [20]. For this 
study, we will choose seasonally adjusted data instead 
of original data. Two classification schemes of the IIP 
are available: indices classified by industry and indices 
classified by use of goods. We adopt the latter classifi- 
cation scheme because we are interested in input-output 
interindustrial relations here, which are measured by cor- 
relations between shipments of final demand goods and 
production of intermediate goods in the IIP data. The 
concept is illustrated in Fig. [T] We emphasize that the 
inner loop of production existing in the economic sys- 
tem may give rise to a nonlinear feedback mechanism 
to complicate the dynamics of the system; outputs are 
reused by the system as inputs for its production activi- 
ties. Table |l] lists the categories [2T] of goods along with 
weights assigned to each of them for computing the av- 
erage IIP. These weights are proportional to value added 
produced in the corresponding categories, and their total 
sum amounts to 10,000. Unfortunately, the resolution 
of the IIP data for the producer goods is quite poor, 
which are just categorized as Mining & Manufacturing 
and Others. 

Figure [2] shows temporal change of the averaged IIP 
data for production, shipments, and inventory during the 
period of January 1988 to June 2009. The ongoing global 
recession is traced back to the subprime mortgage crisis 
in the U.S., which became apparent in 2007. The eco- 
nomic shock has affected Japan without exception, lead- 
ing to a dramatic drop in the production activities of the 
country, as shown in the figure. 

Since some of the entries, such as g = 16 and 17, are 
missing before January 1988, we use the data ^20] for 
the 240 months from January 1988 to December 2007. 
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TABLE I. Classification of goods according to IIP. First, the goods are classified into two categories, "Final Demand" and 
"Producer," and then those two categories are divided into 19 and 2 subcategories, respectively. The central column specifies 
the index g for the subcategories, which has totally 21 values. The number in the parenthesis associated with each species of 
goods shows its weight used to compute the averaged IIP; the weights are normalized so that their total sum is 10,000. 
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FIG. 2. Averaged IIP data Sa for production (thick solid 
line), shipment (thin solid line), and inventory (dotted line) 
as a function of time t. The correlation matrix is calculated 
using the data in the gray shaded area from January 1988 to 
December 2007. 



Further, this chosen period for the study excludes the ab- 
normal behavior of the IIP data due to The Great Reces- 
sion. We denote the IIP data for goods as Sa^g{tj), where 
a = 1, 2, and 3 for production (value added), shipments, 
and inventory, respectively. Similarly, g — 1,2,. ..,21 



denotes the 21 categories of goods, and tj — jAt with 
Ai = 1 month and j = 1,2,...,N{= 240); j = 1 and 
j — N correspond to 1/1988 and 12/2007, respectively. 
The logarithmic growth rate ra,g{tj) is defined as 



Sa,g{tj) 



(1) 



where j runs from 1 to A^' := 
normalized as 



1(= 239). Then, it is 



"a,gV^] 



it. 



r. 



0'a,g 



(2) 



t denotes average over time ti, 
is the standard deviation of r, 



Definition 



over 

2|) ensures that the set ^ := 
■ ■ , Wa,g{tN')} has an average of zero 



where 
and (T, 
time 

{Wa,g{ti),Wa,g{t2), 

and a standard deviation of one 

Figure [3] shows an overview of how the volatility 
ui^ g{tj) of the standardized IIP data behaves on a time- 
goods plane. Unfortunately, the visualization does not 
allow for detecting any correlations involved in the IIP 
data. One may even doubt whether useful information 
on interindustrial relations truly exists in the data. 

To answer the obvious question that would arise here, 
we begin with calculating the equal-time correlation ma- 
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FIG. 3. Bird's-eye view of volatility of standardized IIP data FIG. 4. Temporal variation of strength of fluctuations asso- 
in a panel form, where the index I is defined as I := 21(a — ciated with each eigenmode, where n is an index assigned to 
1) + eigenmodes in descending order of their eigenvalues. 



trix C of {w^ J\ according to 



(3) 



whose diagonal elements are unity by definition of the 
normalized growth rate Wa^g(tj). Since a runs from 1 
to 3 and g (h) runs from 1 to 21, the matrix C has M x M 
(M = 63) components. We denote the eigenvalues and 
the corresponding eigenvectors of the correlation matrix 
as A'-"' and V'-"', respectively: 



(4) 



where the eigenvalues are sorted in descending order of 
their values and the norm of eigenvectors is set to unity. 

On the basis of the eigenvectors V*^"-' thus obtained, 
the normalized growth rate Wa,g{tj) can be decomposed 
into 



n=l 



(5) 



The correlation matrix C is also decomposable in terms 
of the eigenvalues and eigenvectors as 



M 

C = ^ ;^(«)y(")y(n)T^ 



(6) 



The eigenvalues satisfy the following trace constraint: 

M 



J2 A(") = M. 



(7) 



n=l 



By substituting Eq. ^ into Eq. ([s]) and comparing it 
with Eq. ([g]), wc find that 



The eigenvalue of each eigenmode thus represents the 
strength of fluctuations associated with the mode. Fig- 
ure shows the temporal variation of (t) , which is in 




FIG. 5. Probability distribution function p{X) for eigenvalues 
(A) of correlation matrix derived from IIP data, in comparison 
with corresponding result of random matrix theory (RMT) 
represented by solid curve. 



sharp contrast to the results shown in Fig. |3] The trans- 
formation of the base for describing the IIP data reveals 
that very few degrees of freedom actually are responsible 
for the complicated behavior of the IIP. 

Thanks to the RMT, we are able to quantify how many 
eigenmodes should be considered. Probability distribu- 
tion function p{\) for the eigenvalues (A) of the correla- 
tion matrix C is shown in Fig. [5j It is compared with 
the corresponding result [22 of the RMT in the limit of 
infinite dimensions: 




Q V(A+-A)(A-A_) 



for A_ < A < A+, 
otherwise. 



(9) 

(8) where Q :— N' /M ~ 3.79 (> 1) and the upper and lower 
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bounds A± for A are given as 

(1 ± VQr 



A. 



Q 




(10) 



We see that the largest and the second largest eigen- 
values, designated as A^^H- 9.95) and A(2)(~ 3.83), 
are well separated from the eigenvalue distribution pre- 
dicted by the RMT, whereas the third largest eigenvalue 
^(3)(^2± 2.77) is adjacent to the continuum. Therefore, 
only 2 eigenmodes out of a total of 63 are of statistical 
significance according to the RMT. 

Readers may be curious about the present construction 
of a correlation matrix by mixing up data for production, 
shipments, and inventory, because these are very different 
species of data at first glance. Thus far, physicists have 
applied the RMT mainly to analyses of stock data having 
similar characteristics. In this sense our approach is quite 
radical. However, production, shipments, and inventory 
form a trinity in the economic theory for business cycles, 
so that those variables should be treated on an equal 
footing. Using the two dominant eigenmodes, in fact, 
we were successful in proving the existence of intrinsic 
business cycles p/. 

In passing, we note that one may favor the growth rate 
itself defined by 

„ fj- \ Sa^g{tj-^i) — Sa,g{tj) 

'a.gyi'j) ■— Q f-l- \ ^ ' 

for the present analysis over the logarithmic growth 
rate (fTl). If the relative change in Sa,g{tj) is small 



need not distinguish between Eqs. ([ij and (111 numeri- 
cally. To confirm that the results obtained here are insen- 
sitive to the choice of stochastic variables, we repeated 
the same calculation by using Eq. (Ill and found no ap- 



preciable difference between the two calculations for the 
dominant eigenvalues and their associated eigenvectors. 
For instance, the first three largest eigenvalues 9.95, 3.83, 
and 2.77 as shown in Fig. [5] are replaced with 9.96, 3.73, 
and 2.78, respectively. 



III. ROTATIONAL RANDOM SHUFFLING 

There are two major sources of noise in the IIP data. 
One of them, corresponding to thermal noise in physi- 
cal systems, arises from elimination of a large number of 
degrees of freedom from our scope as hidden variables. 
This highlights the stochastic nature of the IIP and has 
a strong influence on autocorrelation of all goods. The 
other source of noise originates from the finite length of 
time series data. Such statistical noise hinders the de- 
tection of correlations among different goods in the IIP 
data. If one could have data of infinite length, statisti- 
cal noise would disappear in the mutual correlations and 
only thermal noise would remain. These two types of 
noise should be distinguished conceptually. The RMT is 



an effective tool for eliminating statistical noise from raw 
data to extract genuine mutual correlations. 

However, the noise reduction method based on the 
RMT heavily depends on the following assumption: 
stochastic variables would be totally independent if cor- 
relations between different variables were switched off. 
Such a null hypothesis simultaneously excludes both au- 
tocorrelations and mutual correlations. In the case of 
daily change in Japanese stock prices that were avail- 
able [17] to us, we found no detectable autocorrelations 
in the corresponding variables; therefore, the RMT func- 
tions ideally. In contrast, the IIP data have significant 
autocorrelations as shown in Fig. [6j where the autocor- 
relation function Ra^g{t) of the normalized growth rate 



Wa.g is defined as 



-^a.g {im ) 



iV' 



N' —rr, 

- V 

TO ^ — ' 



Wa,g{tj)Wa^g{tj+rn)- (12) 



By definition, Ra,g{0) — 1, and if there are no autocor- 
relations, Ra,g{tm) = for TO > 1. We observe that both 
production (a = 1) and shipments {a — 2) have nontriv- 
ial values of autocorrelations at t — 1 month, whereas 
there is no clear evidence for autocorrelations for inven- 
tory (a = 3) in the same time interval; the values aver- 
aged over 21 goods are Ri{l) ~ -0.31, i?2(l) ~ -0.39, 
and Rsi^) — 0.007. Beyond the one-month time lag, 
however, we find no appreciable autocorrelations for any 
of these three categories. 

To formulate the null hypothesis of the RMT using ac- 
tual data, one may shuffle the IIP data completely in the 
time direction. In fact, the eigenvalue distribution of the 
resulting correlation matrix reduces to that of the RMT 
as demonstrated in panel (a) of Fig. [7j where a total of 
10^ samples were generated. Departure from the RMT 
owing to finiteness of the data size is almost negligible 
even for such small-scale data as the IIP. This random- 
ization process inevitably destroys both autocorrelations 
and mutual correlations. From a methodological point 
of view, it is favorable to deal with these two types of 
correlations separately. 

We instead propose to shuffle the data rotationally 
in the time direction, imposing the following periodic 
boundary condition on each of the time series: 



"a,g 



(tj) — >• U'Q,g(tMod(j-T,A"))' 



(13) 



where r € [0, N' — 1] is a (pseudo-) random integer and is 
different for each a and g. This randomization destroys 
only the mutual correlations involved in the data, with 
the autocorrelations left as they are; therefore, it provides 
us with a null hypothesis more appropriate than that of 
the RMT. 

Panel (b) of Fig. [t] shows the result in rotational shuf- 
fling with the same number of samples as that in the 
complete shuffling. We find that the existence of auto- 
correlations alone leads to departure from the RMT. The 
third largest eigenvalue A*^'^-' ~ 2.77 becomes even closer 
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FIG. 6. Autocorrelation functions Ra,g{t) of production {a = 
1), shipments (q = 2), and inventory (q — 3) for each of the 
goods {g = 1,2, ...,21) at t = 1 (a), f = 2 (b), and t = 
3 (c), where the index £ on the horizontal axis is defined in 
the same way as in Fig. [sj The 95% confidence level for no 
autocorrelations is represented by the gray shaded band in 
each panel. 



to the upper limit, A'|_ = 2.47 ± 0.20, of the eigenvalues 
obtained on the basis of the alternative null hypothe- 
sis, where the error is estimated at 95% confidence level. 
This result reinforces neglect of the third eigenmode by 
the RMT. 

Thus, this new method for data shuffling conceptually 
clarifies noise elimination for the correlation matrix, al- 
though the difference in the eigenvalue distribution from 
that of the RMT is practically not very dramatic. In ad- 
dition, we note that the rotational shuffling of the stock 
price data in Japan reproduces the RMT result quite 
well, as is expected from the fact that no appreciable 
autocorrelations are observed there. 
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FIG. 7. Same as Fig.js] but for eigenvalues (A) of the correla- 
tion matrix obtained by shuffiing IIP data completely (a) or 
rotationally (b) in the time direction. It should be noted that 
the autocorrelations involved in the IIP data are preserved in 
rotational shufHing. 



IV. GENUINE CORRELATION MATRIX 

In the current system of IIP data, the above careful 
arguments permit us to adopt 



.(G) 



2 

E 

n=l 



A(")y(")y(«)T _|_ [diagonal terms] (14) 



as a genuine correlation matrix, which consists of just the 
first and second eigenvector components in the spectral 
representation ([6]) of C plus the diagonal terms, thereby 
ensuring that all the diagonal components are 1. We note 
that self-correlations of stochastic variables always exist 
even if they are merely noise. The components of C are 
explicitly written as 



1 



(G) 



for £ = m, 

2 

A(") y/"^ Vl") otherwise. 



(15) 



The eigenvectors V^^-* and V^^-*, associated with A^^-* 
and A(^\ are shown in panels (a) and (b) of Fig. [sj 
respectively. These two eigenvectors have characteris- 
tic features that distinguish them from each other. The 
eigenvector V(^) represents an economic mode in which 
production and shipments of all goods expand (shrink) 
synchronously with decreasing (increasing) inventory of 
producer goods. This corresponds to the market mode 
obtained for the largest eigenvalue in the stock market 
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analyses [HI [TO] , and may be referred to as the "aggregate 
demand" mode according to Keynes' principle of effective 
demand: both shipments and production in all the sec- 
tors are moved jointly by aggregate demand [221 ■ On the 
other hand, the eigenvector V^^) is a mode that appar- 
ently represents dynamics of inventory, i.e., accumulation 
or clearance of inventory, for most goods, including pro- 
ducer goods. We further find positive correlation between 
production enhancement and inventory accumulation for 
most goods. This finding indicates that production has 
a kind of inertia in its response to change of demands. 

Accordingly, we project out raw fluctuations of Wi{tj) 
onto the first and second eigenmodes; that is, only the 
first two terms are retained and the remaining terms are 
regarded as just noise in expansion ([s]): 



'^dt]) = X! ""('^j) V + [noise] 



(16) 



This process extracts statistically meaningful informa- 
tion on mutual correlations among wiitj) as has been 
already discussed in Sees. |ll] and |III[ Collaboration of 
these two modes results in inherent business cycles with 
periods of 40 and 60 months throughout the economy. 
The cycles are accounted for by time lags in information 
flow between demand of goods by consumers and decision 
making of firms on production [1] ; inventory fills this in- 
formation gap. If we singled out the most dominant mode 
alone in Eq. ([5| , all wi {tj ) would oscillate without phase 
difference. As will be shown later, each goods possesses 



its own characteristics in the phase relations among pro- 
duction, shipments, and inventory. 

Temporal change of the two principal factors ai(t) and 
02 (i) is plotted in panel (a) of Fig.joj Since the functional 
behavior of these variables is very noisy, we take their 
simple moving average defined as 



anitj) := 



2C + 1 



(17) 



where ^ is a characteristic time scale for smoothing. This 
process eliminates "thermal" noise present in the origi- 
nal data. Actually, the moving average was taken with 
^ = 6; the results for ai{t) and 02 (t) are shown in panel 
(b) of Fig. [9] We see that the moving-average operation 
significantly reduces the level of noise present in a„(i). It 
is noteworthy that Fig. [9] indicates the existence of some 
mechanical relationship between ai(t) and 02 (t). This 
finding is ascertained more quantitatively from Fig. |10[ 
in which the correlation coefficient between ai{t) and 
a2{t - r), 

a2it-T)) 

^ir) = ; ^ , (18) 

is plotted as a function of time lag t. A correlation as 
large as 0.7 is detected between the two dominant modes 
around r = 10 months. Detailed study of the underlying 
dynamics in the economic system is in progress and will 
be reported elsewhere. 



V. FLUCTUATION-DISSIPATION THEORY 

The fiuctuation-dissipation theory plays a central role 
in nonequilibrium statistical mechanics because this the- 
ory establishes a general relation between fluctuation 
properties of a physical system in equilibrium and re- 
sponse properties of the system to small external pertur- 
bations. We assume that the theory still is applicable to 
the economic system under study here. This assumption 
provides us with a framework to derive input-output in- 
terindustrial relations in the system. Its validity in view 
of how the system responded to the recent economic crisis 
will be discussed later . 

Let us denote our variable Wa.g as Wi, whose index 
£ :— 21{a — 1) + g runs from 1 to 63. We assume that wi 
obeys dynamics governed by a Hamiltonian H{{w}, {a;}), 
where {x} is a set of "hidden" variables in the system, 
which encompass all variables in the current economics. 
These numerous variables interact with each other in a 
nonlinear chaotic way and hence, the temporal change 
of We appears stochastic in the same way as that of a 
Brownian particle. The existence of underlying dynamics 
in the IIP 1^, as demonstrated in Figs.[9]and 10 strongly 
supports this idea borrowed from mechanics of motion. 
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FIG. 10. Correlation coefficient C^y a^(r) between a\{t) and 
02 (t — t) as a function of time lag r, calculated during the 
normal period (from January 1988 to December 2007) with 
^ = 6 for the moving-average operation. A comparison of this 
result with that obtained using unsmoothed data with ^ = 
is also shown in the figure. 



Actually, however, the economy of a nation is quite 
open now; therefore, it could potentially be subjected to 
perturbations such as disasters, political issues, and trade 
issues. We thus add external forces tg{t) to the system; 
then, the total Hamiltonian T-L becomes 



M 



n = H{{w},{x}) -Y,e,{t)w,. 



(19) 



This extra term represents external perturbations to the 
equation of motion for wf. 



dt 
dwi 
~dF 



dH 
dwi 
dH 
dpi' 



(20) 
(21) 



where pi is the momentum conjugate to wg. Therefore, 
ei{t) directly affects wg at time i, the effect of which 
then extends to other w's through direct and indirect 
interactions among them. 

For simplicity, let us assume that {e} is constant in 
time. Thus, the perturbation set {e} induces a static 
shift of the equilibrium positions of the variables {w}, 
which otherwise move stochastically around the origin 
wi — 0. If the perturbation is weak, the shift {wi) thus 
induced can be expressed by the following linear response 
relation: 



M 



(22) 



where the ensemble average denoted by (•) has replaced 
the time average. The coefficients {%} are the result of 
the interactions, and they are called "magnetic suscepti- 
bility" while describing the physics of magnetic materials. 

Once such a set of susceptibilities is available, we can 
quantify the response of the economic system to external 
perturbations. For instance, suppose that the govern- 
ment adopts an economic policy to increase the shipment 
of one of the final demand goods by {wm) with a stimulus 
tm- The resulting changes in production, shipments, and 
inventory of goods are given as 



(Wl) = Xlr 



{wm) = XMr, 



(23) 



Since Cm is not an observable quantity, it should be ap- 
propriate to eliminate Cm appearing in Eq. (23) and ex- 
press ripple effects on the economy in terms of {wm) as 



^1 



Xlr. 



(24) 



/ \ XMm I 

(wm) = {w„ 
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In Sec. VI we demonstrate that Eq. (24 1 can provide 



quantitative information on input-output interindustrial 
relations in Japan. Further, one may make reverse use of 
the hnear response relation ( 22 ) to distinguish and detect 



external perturbation from observed economic changes in 
{w}; this is discussed in Sec. 



VII 



Now, the remaining problem is how to calculate {x}- 
To this end, we invoke the concept of the fluctuation- 
dissipation (FD) theorem in statistical physics. If we 
assume that the stochastic process of {w} is characterized 
by Gibbs' ensemble, then the probability density function 
(PDF) Pi{w}, {e}) for {w} is given as 



P(M,{e}) (X exp[-/3H(W,{e})], 



(25) 



where the hidden variables {x} have been integrated out 
and f3 is the inverse "temperature" of the economic sys- 
tem. For a weak perturbation, Eq. ( 25 ) is expanded to 
the first order of {e} as 



Pi{w}, {e}) P{{w}) + E ' 



where 



P{{w}) cx exp[-/3H{{w})] 



(26) 



(27) 



is the PDF in the absence of {e}. Equation (26 1 enables 
us to calculate the induced change in W£ by {e| as 



(wi) 



P{{w},{e})wid{w} 



^ E ^™ / P{{'w})we'Wmd{w} (^28) 

m=l 
M 



where denotes the ensemble average without pertur- 
bation. Comparison of Eqs. (22) and (28) gives one of 



the outcomes of the FD theorem: 



(29) 



where C'^*^-' denotes a correlation matrix in the absence 
of external perturbations. 

Making use of Eq. ( 29 ) , we can rewrite the relation 



( 24 ) of economic ripple effects caused by the increase in 



shipments of final demand goods as 



(30) 



Wm) = Cfllni^ra)- 



This is just one example of possible interindustrial re- 
lations derived from the present formulation. What we 
should emphasize here is that Eq. ( 30 ) has a rather gen- 



eral form in the framework of linear response. 



For instance, we do not need to determine the temper- 
ature of the economic system for /3. Even the assump- 
tion of Gibbs' ensemble, e.g., as given in Eqs. (25) and 



(27), may be too restrictive, because the assumption (26) 



about the PDF is sufficient to derive Eq. ( 30 ) . We also re- 
call Onsager's regression hypothesis [24ll25j . on which the 
fluctuation-dissipation theorem relies. Once one accepts 
the hypothesis, one can readily derive Eq. ( [M] ). Accord- 
ing to him, the response of a system in equilibrium to an 
external field shares an identical law with its response to 
a spontaneous fluctuation. In other words, the regression 
of spontaneous fluctuations at equilibrium takes place in 
the same way as the relaxation of non-equilibrium dis- 
turbances does. Let us suppose that the non-equilibrium 
disturbances (wi) and are linearly related through 



(Wi) = K{Wm)- 

Accordingly, the spontaneous fluctuations Wi and Wn 
isfy the same relation as Eq. (31): 



Wi 



KW„ 



(31) 



sat- 



(32) 



The ensemble average of Eq. ([32j) multiplied by w„i on 
both hand sides determines the proportionality coeffl- 
cient K as 



(0) 



(33) 



We thus see that Eq. ( 30 ) is directly derivable from On- 
sager's hypothesis. 

We note that the correlation matrix appearing in 
Eqs. (29) and (30) should be measured for a system not 



subject to any perturbations. However, t he g enuine cor- 



(G) 

relation matrix C) determined by Eq. (15) is possibly 



contaminated with various kinds of external economic 
shocks. While such forces may easily affect the stochas- 
tic motion of each w, it is legitimate to assume that their 
influence on the correlations among w's are much weaker; 
otherwise, the external factors would have to work coher- 
ently to change "springs" connecting pairs of w's. This 



consideration justifies the replacement of ci'fi in Eq. (|30|) 



with C, 



(G) 



VI. INTERINDUSTRIAL RELATIONS 

We are now in a position to quantitatively estimate 
the strength of the interindustrial relations by making 
use of the genuine correlation matrix through Eq. ( 30 ) . 
In particular, we focus on ripple effects on production 
of intermediate goods that are triggered by applying an 
external stimulus to consumption of final demand goods: 



{Wl,2l) = C- 



_ ^(G) 



l,21;2,g 



(34) 



(35) 
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FIG. 11. Input-output interindustrial relations based on the 
genuine and original correlation matrices. The vertical axis 
in panel (a) indicates the extent to which the logarithmic 
growth rate of production of intermediate goods for Mining 
& Manufacturing {g = 20) is accordingly elevated when the 
logarithmic growth rate of shipments of each of the final de- 
mand goods specified on the horizontal axis is increased by 
one unit. Panel (b) is the same as panel (a), except that 
it shows the relationship between production of intermediate 
goods for Others (g — 21) and shipments of each of the final 
demand goods. 



with g = 1, 2, • • • , 19. 

The results for production of intermediate goods for 
Mining & Manufacturing {g — 20) are shown in the up- 
per panel of Fig. [TT] in which those obtained with the 
original correlation matrix are also added for comparison. 
This figure shows an increase in the logarithmic growth 
rate of production of intermediate goods that is predicted 
from unit increment of the logarithmic growth rate of 
shipments of each of the final demand goods. As ex- 
pected, increase in shipments of final demand goods with 
large weights, as represented by Manufacturing Equip- 
ment {g — 1), Construction (g = 9), Motor Vehicles 
(g = 15), House Work {g = 16), and Food & Beverage 
{g = 19), certainly causes large ripple effects on the pro- 
duction of intermediate goods. If the original correlation 
matrix is replaced with the genuine one, then the rela- 
tive importance of species of final goods is interchanged 
between Construction and Motor Vehicles. This is under- 
standable because sales of cars are sometimes promoted 
just for inventory adjustment, having no effect on the 



growth of production of intermediate goods. We also note 
that the original correlation matrix significantly under- 
estimates the effects of Furniture & Furnishing {g = 13) 
and Nondurable Consumer Goods {g = 17,18,19). It 
is noteworthy that Furniture & Furnishing and Clothing 
& Footwear, having much smaller weights than the ma- 
jor final demand goods, have comparable contributions; 
some feedback mechanism must be working through the 
inner loop in the economic system. 

The lower panel in Fig. [Tl] shows the corresponding 
results for production of intermediate goods for Others 
{g = 21), whose weight is one order of magnitude smaller 
than that of Mining & Manufacturing. The important 
species of final demand goods are common in both cat- 
egories of intermediate goods. In contrast, the original 
correlation matrix significantly underestimates the effects 
of different final demand goods such as those given by 
g = 3 to g — 7. 

Presence of a correlation between two stochastic vari- 
ables A and B does not indicate the existence of a me- 
chanical connection between them. Actually, correlating 
A and B might be driven by a third variable C; then, 
there would be no causality relationship between A and 
B. To address this question, we provide detailed informa- 
tion on phase relations in the business cycles identified 
in the previous study [T . Table |ll] lists phases of the 
cyclic motion of production, shipments, and inventory at 
T = 60 and 40 for each of the goods. We can see that 
shipments of final demand goods are ahead of or almost 
in phase with production of intermediate goods; the res- 
olution limit (one month) is 6° and 9° for T = 60 and 
40, respectively. Electricity {g = 2) and Communication 
& Broadcasting {g = 3) are exceptions to this obser- 
vation. The wave of production arrives first and then, 
that of shipments follows for Electricity; the production 
activity for Communication & Broadcasting behaves sig- 
nificantly out of phase with that averaged over goods. 
Since we have adopted a static approximation for the in- 
terindustrial relations, it may be more appropriate to av- 
erage the phase relations over frequency. The results are 
shown in Fig. [12] and Table |111| The frequency-averaged 
phase relations in the cyclic behavior of the economic 
fiuctuations thus support our postulate that production 
of intermediate goods is driven by increasing shipments 
of final demand goods with a few exceptions. 



VII. EXTERNAL STIMULI 

Finally, we try to identify the presence of external stim- 
uli hidden in real data by inversely using the linear re- 
sponse relationship (22). The recent global economic cri- 



sis certainly has delivered an extremely large shock to the 
economic system of Japan, as is clearly shown in Fig. [2] 
In our previous paper [1 , however, we demonstrated that 
the crisis has simply increased the level of fluctuations as- 
sociated with the dominant modes that were determined 
from the data during the normal time, instead of destroy- 
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TABLE II. Phases of periodic oscillations with T = 60 and 
40 of production (P), shipments (S), and inventory (I) for all 
goods, in the unit of degrees ranging from —180 to 180. They 
are measured relative to the production of = 20. 
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19 
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21 
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Average 


-0.6 


25.2 


-66.6 


-0.6 


16.9 


-123.7 



TABLE III. Frequency-averaged phases of periodic motion of 
production (P), shipments (S), and inventory (I) for final de- 
mand and producer goods in the unit of degrees. The results 
for final demand goods were obtained by averaging over goods 
excluding g = 2 and 3; the numbers in parentheses are those 
obtained with all final demand goods. 



Final demand goods 



Producer goods 
g = 20 5 = 21 



P 
S 
I 



-0.69 (-2.35) 

2.99 (0.20) 
-28.7 (-29.5) 




5.84 

-37.4 



0.99 
5.15 

-37.5 



ing the industrial structure itself; this is also manifested 
here, as shown in Fig. [9] And the information on collec- 
tive movement of the IIP that we could extract from the 
dominant modes remains intact even in such an abnor- 
mal situation. This result thus conforms to the idea of 
Onsager's regression hypothesis, indicating the validity 
of the fluctuation-dissipation theory even in an economic 
system that is supposed to be far away from equilibrium. 



Since approximation ( 15 ) has been adopted for the cor- 
relation matrix, we consider only two independent exter- 
nal fields {771, 772} that are coupled to the normal coor- 
dinates {ai, 02} associated with the two dominant eigen- 




FIG. 12. Frequency-averaged phases of production, ship- 
ments, and inventory for each goods, measured relatively to 
production of 3 = 20 in the unit of degrees. 



modes V^^^ and V'-^\ respectively. The total Hamilto- 



nian ( 19 ) is therefore simplified to 

H = (fll, 02, {x}) - yyiai - 77202. (36) 



The reduced external fields {77} in Eq. ( 36 1 are derived 



from the original ones in Eq. ( 19 1 through 

M 



77n 



(") 



(37) 



Then, Eq. (22) is projected onto the two-dimensional re- 



duced state space as 



{ai) 

(«2> 



where 



Xii X12 
X21 X22 



M 



(a„) = 5] (u;,) 



(38) 



(39) 



e=i 

and the reduced susceptibilities {x} are defined as 

M M 

1=1 j=i 



(40) 



The relative values of {x} with reference to xn are cal- 
culated from C*-*^-* as 



Xii X12 
X21 X22 



/3 



1 1.30x10-3-, 
1.30 X 10-3 0.433 I ■ ^ > 



This results shows that the two eigenmodes are almost 
decoupled from each other, which is understandable from 
the orthogonality (|8| of the normal coordinates. 



One can obtain {r/} using the inverse of Eq. ( 38 ) along 
with Eqs. (39 1 and (40|, although it is not so straightfor- 
ward. We first recall that (wf ) in the right-hand side of 
Eq. ( |39[ ) is the deviation of wi from the equilibrium value 
induced by external perturbation, and not fluctuations 
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Figure ITS) for which we arbitrarily set /3 = 1 in 



. 1.0 I I I I I I I I I I I I I I I I I I I I I L 

1988 1990 1992 1994 199G 1998 2000 2002 2004 2006 2008 2010 
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FIG. 13. External stimuli 771 and 772, derived from IIP data 
through the linear response relation ( 38 1 , shown as a func- 



tion of time in panels (a) and (b), respectively; the system is 
assumed to respond instantly to the applied external fields. 
The solid curves depict results obtained only with the terms 
of fc = l,2,4,and6 in Eq. (42l, and the dotted curves depict 
those calculated with the terms of fc ^ 9. The gray shaded 
area is the same as depicted in Figs. [2] and [5] 



of We directly observed in the real data. We then iden- 
tify (wi) as residuals obtained by subtracting the long- 
period components arising from the inherent business cy- 
cles from moving average fluctuations of the IIP. 

To extract (we), we first define the Fourier transform 
of the coefficients a„(ij) as 



Eq. (41 ), shows the external fields ryi and 772 thus derived 
from (u!^). Two computational schemes were adopted to 
evaluate the long-period components in the IIP data, and 
no appreciable difference was observed between the two 
results. Here, the economic system was assumed to re- 
spond instantaneously to the applied external fields with- 
out any time delay. Referring to Fig. [2] we clearly con- 
firm that such a large external shock as manifested in 771 
causes the drastic drop in industrial activities in Japan. 
We also see that another large shock in 772, which leads 
to reduction in inventory, immediately accompanies the 
first shock. In contrast, the maximum fluctuation levels 
of 771 and 772 are 0.1 and 0.2, respectively, in the normal 
period (before the end of 2007). 



VIII. CONCLUSION 

This paper has described our attempt to utilize the 
fluctuation-dissipation theory for elucidating the nature 
of input-output correlations in the Japanese industry on 
the basis of IIP data. We were able to quantitatively esti- 
mate the strength of correlations between goods by using 
the genuine correlation matrix obtained in this study. We 
were also successful in extracting external stimuli over 
the last two decades. The noise reduction along with 
the RMT enabled us to detect economic signals hidden 
behind the complicated dynamics of the IIP. The strong 
coincidence between the sudden change in IIP data and 
the external shocks described here may prove that the 
present method is capable of predicting the input-output 
interindustrial relationship with a much higher time res- 
olution than the annual resolution. We thus expect the 
results of this study to provide a new methodology for 
gaining deeper understanding of complex economic phe- 
nomena at a macroscopic level. 
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JV'-l 



k—1 



an{u)k) e 



'iuiktj 



(42) 



with the Fourier frequency ujk — 2nk/{N'At) and hence 
LUktj = 2nkj /N' . The relevant long-period component 

Om^^tj) is obtained by limiting the sum over k only to 
k = I {T = 240), fc = 2 (T = 120), /fc = 4 (T = 60), 
and fc = 6 (T = 40) or by summing all of the terms with 
periods larger than 2 years {k ^ 9) in Eq. (42). The 
formula for {wi) is finally expanded as 
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